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Abstract 

We develop method that allows to derive reductions and solutions to hyper- 
bolic systems of partial differential equations. The method is based on using 
functions that are constant in the direction of characteristics of the system. 
These functions generalize well-known Riemann invariants. As applications 
we consider the gas dynamics system and ideal magnetohydrodynamics equa- 
tions. In special cases we find solutions of these equations depending on some 
arbitrary functions. 

1 Introduction 

One of the first methods for finding solutions to nonlinear partial differential 
equation 

P 1 y 1 ^ 1 Uy , Uxx 1 Uxy, Uyy ) — (1) 

was proposed by Monge and was then further improved by Ampere. The 
method was described in detail in the classical books of Goursate ^ and 
Forsyth [2]. To apply this method, one must find an equation of the first 
order 

f{x,y,u,Ux,Uy) = c, ceR (2) 

such that every solution of ^ satisfies the equation (pQ) for arbitrary c. In this 
case the function / is called a first integral of the equation To find first 
integrals, we need to look for functions which are constant in the direction of 
characteristics of the equation . If there are two first integrals /i and /2 for 
given family of characteristics, then integration of the equation reduces to 
solving of the first order equation 

G(/i,/2)=0, 

with G an arbitrary function. 

In 1870, G. Darboux ^ announced a generalization of the Monge- Ampere 
method. He proposed to seek an additional partial differential equation of 
second order (or higher) 

dip^ iV I'^i'^x t'^yt'^xxi'Uxyt'^y'y) — C (3) 
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such that the system of equations and Q is in involution for all c. The 
function g turns out to be constant along a family of characteristics of the 
equation In this case, the function g is called a characteristic invariant 
of When partial differential equation (pQ) has a sufficient number of char- 
acteristic invariants, it can be reduced to a system of ordinary differential 
equations. The detailed description of the Darboux method is also given in 
the above mentioned books [J and [Ij. 

Although equations that are integrable by Darboux's method arise rarely, 
they are of much interest. E. Vessiot |5j classified all equations of type 

Uxy = w{x,y,u,Ux,Uy) 

integrable by this method and found a general solution for every obtained 
equation. Recently, there was a renewed interest in the method of Darboux 
that was studied in 0, d, [H], 0, CHI and [TT] . 

In this paper we consider systems of partial differential equations in two 
and n independent variables. In section El we introduce an operator of differ- 
entiation in the direction of characteristics of the system and corresponding 
invariants of characteristics of order k. We prove that if a function h defined 
on the k-th order jet space J^^^ is constant along a vector field on the solu- 
tions of the system of partial differential equations, then this function is an 
invariant of characteristics. In section |21 we describe the Darboux method 
for systems of partial differential equations in two independent variables and 
give its applications to the gas dynamics system and magnetohydrodynamics 
equations. We use invariants of characteristics to reduce these systems and 
find solutions depending on some arbitrary functions. 



2 Invariants of characteristics 

Let us begin with a system of first order partial differential equations in two 
independent and m dependent variables 

ut + F{t,x,u,Ux) = 0, (4) 

where u = (u^, . . . , u™), ut = (n^,...,u™), = {u\, . . . ^u^)^ and F = 

We denote by Dt and the total derivatives with respect to t and x. 
Consider differential operator 

Dt + \Dx, (5) 

whose coefficient A can depend on t, x^ u, and Ux- The operator @ is an 
operator of differentiation in the direction of characteristics of the system , 
if the coefficient A satisfies to the equation 

d{F) _ 

where jt/ — r — 1 mV Jacobi matrix and E is the identity 

0{Ux) 0[Un^,...,Un^ ) 

matrix. 
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Suppose tt^ is the partial derivative of order k of the function with 
respect to x, then Uk = (uj,, stands for the vector composed of these 

derivatives. Let L be an operator of the differentiation in the direction of 
characteristics of the system (jlj). According to .11^, a function h{t, x,u, . . . , Uk) 
defined on the /c-th order jet space j'^^'^ is called an invariant of characteristics 
of order k of the system Q corresponding to the operator L, if /i is a solution 
of the equation 

m\[s\ = ^- (7) 

Here [S] means the system @ and its differential consequences with respect 
to X. When the system ^ has the Riemann invariants, they are zero order 
invariants of characteristics. 

Some systems have invariants of characteristics of arbitrary order. For 
example, consider one-dimensional system of gas dynamics equations |12j : 

Ut + uu^ +Px/P = 0, 

Pt + {pu)x = 0, (8) 

St + USx = 0, 

where u, p, p, and s are the velocity, the density, the pressure, and the entropy. 
The equation of state is given by the function p = p{p, s). The operators of 
the differentiation in the direction of characteristics of the system Q are 

Li = Dt + uD,^, L2 = Dt + {u + c)D,^, Ls = Dt + {u-c)D,^, (9) 
/dp 

with ^ ~ \ / ^ speed of sound. Obviously, the entropy s is an invariant 
of characteristics corresponding to the operator Li. It is easy to check that 

the operator —D^ commutes with Li by virtue of the second equation of the 
P 

system (jSJ. This implies that the recurrent formula 

In+i = -Dx{In), n = 0, 1,... 
P 

gives the invariants of characteristics corresponding to the operator Li. We 
will show in section that invariants of characteristics corresponding to the 
operators L2 and -L3 exist only for the special equations of state. 

It can be proved that if hi and /12 are invariants of characteristics 
of the system corresponding to the operator L, then both an arbitrary 

D h 

function /(/ii, /12) and h = invariants of characteristics. 

Lemma 1. Let L he an operator of the form ©. Suppose that a function 
h{t,x,u,ui, . . . ,Un), with n > 1, satisfies 0, then L is the operator of the 
differentiation in the direction of characteristics of the system @ and h is an 
invariant of characteristics corresponding to the operator L. 

Proof. According to condition of the theorem, /i is a solution of the equation 

Dth + XD,h\is] = 0. (10) 
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Note that 

m 

Dxh ^^K+iKi^, 
1=1 

where the symbol ~ means that the difference between left and right-hand 
sides contains no derivatives of order greater than n. It is easy to see that the 
formula 



l<i,j<m 

is correct because of the system Q). From the equation H10() we have 



m 



l<j,j<m i=l 

This yields m equations 

m 

i=i 

with (5* the Kronecker symbol. Rewriting the above equation in the matrix 
form 

where (/i^i , . . . , /i^™)* is the transposed vector, we conclude that A is a solution 
of the equation □ 

We now consider the system of first order partial differential equations in 
n + 1 independent and m dependent variables 

+ F(t,X,M,Mi.j, . . . = 0, (11) 

with X = (xi, . . . , u = {v} , . . . ,vJ^), = {ul.., . . . ,u^.), and F = 
{F\...,F^). 

Let us denote by the set of k-th order partial derivatives of the functions 
v},...,u^ with respect to xi,...,Xn- We say that a function 
h{t,x,u,ui, . . . ,Uk) defined on the k-th order jet space is an invariant 
of characteristics of the system (|11|) if h satisfies to the equation 

^ djF) 



det i^EDth + Y dfA^'^'V '[^"^ " ^^^^ 



Here Dt and D^^ are the total derivatives with respect to t and Xj, [Sn] 
means the system (|12j) and its differential consequences with respect to Xi 



d(F) d(F^ , . . . , F"^) 
(i = l,...,n), -— r = ' ^ is the Jacobi matrix and E is an 

identity matrix. 
An operator 

L = Dt + XiDx, + --- + XnDx^, 

where Aj can depend on t,x,u,ui, is called an operator of differentiation in 
the direction of the vector field i; = (1, Ai, . . . , Xn)- 
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Theorem 1. Suppose that there is an operator L of differentiation in the di- 
rection of the vector fieldv = (1, Ai, . . . , A„,) and a function h{t, x,u,ui, . . . , Uk), 
with k > 1, such that 

L{h)\[Sn]=0, 

then h is an invariant of characteristics of the system . 
Proof. Since h satisfies 

n 

Dth + Y,^iD.M[Sn] = 0, (13) 
1=1 

the coefficients of n+ 1-th derivatives on the left-hand side of H13|) must vanish. 
To find these coefficients, we write D^^h up to k-th derivatives: 

m m 
j=l \a\=k j=l \a\=k 

Here Ua denotes the derivative — q,^ of order k = \oi\ = ai + ■ ■ ■ + 

is the derivative — — , and the symbol ~ means that 

OX^ ^■■■OX- ' ■■■OXji'^ 

the difference between the left- and right-hand sides contains no the k + 1-th 
derivatives. This yields 

n m 

XiD^.h + ■■■ + XnD^^h -Y^YY^ K+uKi- 

i=l j=l \a\=k 

On the other hand, we have 

m m n m 

msni - -EE ^''inKi - -EE (EE<+i.^O^- 

j=l \a\=k j=l \a\=k i=l s=l 

The above calculations lead to 

n m n m 

Dth + YKD^MiSrA^-Y. E [E(E<-M.^4|-^^<+i.)]^^.=0- 

i=l j=l \a\=k j=l s=l 

It is convenient to represent the last relation in a matrix form 

n 

E E ^"+1.^"' V = 0, (14) 

i=l \a\=k 

With Ua = • -,«„), = (/ini • ■,K^), and A^' = 'q(J[^ -XiE. 

We need to prove that /i is a solution of the equation ((T^ which is equiv- 
alent to the following: 

n 

det(E^"'^x.^) =0- (15) 

i=l 
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For this purpose, it is enough to show that the hnear homogeneous system 

n 

[Y,A^^D,My = (16) 

i=l 

has a nontrivial solution r. This solution is expressed in the form 

r=Y,{DhrK^, (17) 

\a\=k 

where {Dh)°' = {D^^h)'^^ ■ ■ ■ (D^^/i)"". Indeed, substituting ^ in the left- 
hand side of ()13() leads to 

n n 
1=1 \a\=k 1=1 \a\=k 

Note that, the expressions including Ua+u in (|14j) coincide with ones including 
{Dh)°'+^i in ((IHI). Since the left-hand side of (UU is zero then ((THl) is equal to 
zero as well. Hence, (|15|) is valid. □ 

When the conditions of the theorem are satisfied, we say that a function 
h defined on the A;-th order jet space J^*^^ is constant along a vector field 
V = (1, Ai, . . . , An) on the solutions of the system (fTT|). 



3 The Darboux method and its apphca- 
tions 

In this section, we will describe the Darboux method for systems of partial 
differential equations and give relevant examples. The detailed description of 
applications of this approach to second order partial differential equations in 
two independent variables is given in pp. 

The Darboux method is based on using the invariants of characteristics. 
Let us consider the system (jlj and assume that the corresponding equation 
© has m distinctive real roots Ai, . . . , A^. If there are two functionally inde- 
pendent invariants of characteristics Jj for all Aj, then we can constitute a 
system of ordinary differential equations in the independent variable x 

/l(/l,Jl) = 0, fm{Im,Jm)=0, (19) 

where fi,...,fm are arbitrary functions. It is necessary to look for the general 
solution of the systems Q and (fT9|) . If we can find a general solution of (fT9|) , 
then substituting this solution into (j^ leads to a system of ordinary differential 
equations in the independent variable t. It is enough to solve the last system 
in order to find the general solution of (jl)). 

Remark. When the equation Q has a root of multiplicity k it is desirable to 
obtain k + 1 invariants Ji, . . . , Jk+i corresponding to this root. In this case, 
the system (fT^ includes equations 

/l(J2, Jl) =0,...,/fc(Jfc+i,Ji) =0. 

Such example arises naturally in the equations of magnetohydro dynamics. 
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As example, consider the system of gas dynamics equations in two inde- 
pendent variables t and x: 

Ut + UUx +Px/p = 0, 

pt + (pu)^ = 0, (20) 

Pt + UPx + PC^Ux = 0. 

where p is density, u is velocity, p is pressure, and c{p,p) is sound speed. 

One can easily deduce zero order invariants of characteristics of the system 
(PU)) corresponding to the operator L2 (or L3) given by ©• To do so, following 
[TT] , one must seek all solutions of the equation 

Dth + iu + c)D,h\[G]=0, (21) 

where [G] stands for system H20() and its differential consequences with respect 
to x; the function h can depend on t,x,u,p,p. Obviously, the left-hand side 
of (PT|) is a polynomial of the first degree in Ux, Px, and px- Collecting similar 
terms of these variables leads to the following equations 

hp = 0, hu = pchp, ht + {u + c)hx = 0. (22) 

It follows from the first and second equations of H22|) that a nonconstant 
solution h exist only if 

c = 9{p)/p, (23) 

with g an arbitrary function of p. As a consequence of the third equation, h 
is independent of t and x. According to the second equation of (|22|) . h is an 
arbitrary function of 

dp 

gip)' 

Similarly, it is possible to check that the Riemann invariant 



r = u + 



I~ =u 



dp 

g{p) 



corresponds to the operator L3. 

We now use the invariant /~ to derive solutions of the system H2U|) . Setting 

f ^P 

1^=0 and introducing a new function F = , we get the following 



9{p)' 



representation 

u = F{p). 

In this case, the system reduces to 

pt + {Fp)x = 0' Pt + Fpx + j^ = 0. (24) 
The previous system admits the Riemann invariants 

p, r = J {F'fdp+^. 
One then can rewrite the system ()24p in the form 

Pt + (^F{p) - ^^1^) Px = 0, rt + F{p)r, = 0, 
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where G{p) = J (^F' {p))^ dp. Using the hodograph transformation leads to 
the hnear system 



F{p)tp = 0, {F{p)F'{p) - Gip) + r)tr - F'{p)xr = 0. 



The general solution of this system is 



t = P + R'F\ 



R' {FF' + l/p) -R + j FP'dp, 



(25) 



where P = P{p) and R = R{r) are arbitrary functions. As a consequence of 
(|25() . we obtain a solution of (|2U() in the implicit form. 

It was shown in that the first order invariants corresponding to the 
operators L2 and L3 exist only if the speed sound is given by 



c=(a + 6p)(2/3)/p, 
The corresponding equation of state is 



a,b e R. 



pip, s) 



a + 



3/0 



b{A{s)p-l] 



the Riemann invariants are 

/2 = 6u + 3(a + I3 = bu - 3{a + bp) 

and the first order invariants have the form 



1/3 



J2 



p{a + bp) 



1/3 



bt 



J3 



u 



p{a + bp)^/^ 
t{a + bp)'^/'^ -Px 



J' 



The invariants corresponding to the operator Li = Dt + uDx (mentioned in 
the section 121) are 



h = b/p - 3(a + bp) 



'1/3 



Ji = 1/pDMi)- 



The gas dynamics system is conveniently written in terms of the Riemann 
invariants 



(^i)t 
ih)t 



h + h 



2b 

hUl-il 



(I- 



l)x, 



2/2 M 



(I: 



3)t 



36b 

hill - li) - 2hM 



36b 



ih)x, 

ih)., 



(26) 



with M = Ii{l2 — I3) + 18. The first order invariants Ji,J2, and J3 can be 
represented as 



M 



h-h 



ih)., 



■h = t 



18b 



Mil 



J 3 = t 



18b 



2)x 



Mil 



■i)x 



(27) 



We now apply the Darboux method to reduce the system (|26j) to some 
ordinary differential equations. The corresponding system ()19() is equivalent 
to 

Ji = Fiih), J2 = F2ih), J3 = Fsih), (28) 
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where Fi, Fi, and F^ are arbitrary functions. Prom (|26|) . (|27|) and (|28j) we get 
two systems of ordinary differential equations: 

... _ Fl(Il)(/2-/3) 

186 

^^^^^ = M(t-F2(l2))' 

M(t-F3(/3))' 

^_Fi(7i)(7|^ 
2M6 
-/|) -2I2M 

- 2M(t-F2(/2)) ' 
/i(/|-/|)-2/3M 



2M(t-F3a3)) ■ 

Introducing new functions 

^(/i) = y -^^dh, G,{Ii) = I F,{h)dh, i = 2,3, 
one may write the equations H29() in the fohowing way 

[hx - ^{h)], = i^, [th - G,{h)], = 1^, [th - G,{h)], = ^. 
Hence, the system (|29() has the first integrals 

tl2-G2{l2)-bx + ^{h)=C2{t), 

th-Gs{h)-bx + ^{h) = C3{t), 

with C2{t) and C3(t) arbitrary functions. Differentiating the previous relations 
with respect to t and using the system (|5n|l . we deduce that the functions C2{t) 
and C3(t) are constants. 

Therefore, the system (|26p can be reduced to a couple of differential equa- 
tions 

bhih-h) 



(I 



l)x 



^'(/l)(/l(l2-/3) + 18)^ 



,.,,2 (31) 
2vI/'(/i)(/i(/2-/3) + 18)' 

where I2 and Is must be expressed from relations 

th - G2{l2) -bx + ^(/i) = 0, th - Gsih) -bx + ^{h) = 0. 

It is possible to find solutions of gas dynamics equations by integrating the 
equations (pTT|) with partial functions G2-,G^, and 

It is interesting to note that there are the second order invariants of char- 
acteristics 

^'"^ " 5 ^ {p'l^u.±Pxf ' 

, _ Gp^{pu^G±G,) 
(2) - 
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with G = , corresponding to the operators = Dt + (u i: c)Dx, 

p^Ux ± px 

when the speed sound is given by one of the fohowing formulas 

c = p'^/^/p, c = p^/p. 
We now consider the one-dimensional magnetohydrodynamics equations 

m 

Pt + {pu)x = 0, 

ut + uux H \ = 0, 

p 47rp 

vt + uvx = 0, (32) 

wt + uwx = 0, 
{B2)t + {uB2)x = 0, 
{B^)t + {uB^)x = 0, 

St + USx = 0. 

Here p is the density, p is the pressure, and s is the entropy; {u,v,w) and 
{Bi, B2, B^) denote the velocity and the magnetic fields, respectively. We 
assume that Bi = and p is a function of p and s. In this case, there are the 
following invariants of characteristics 

T T T J B2 B3 Sx 

h = s, 12 = V, h = w, 14 = — , I5 = — , J = 

P P P 

corresponding to the operator Li = Dt + uDx- 

Using these invariants we will reduce the system (|32|) to the second order 
equation for the entropy. According to the general scheme of the Darboux 
method, we write 

v = V{s), w = W{s), Sx/p = fi{s), B2/p = f2{s), Bs/p = f3{s), 

where V, W, /i, /2, and /s are arbitrary functions. These relations are equiva- 
lent to the following representation 

v = V{s), w = W{s), p = s,n(s), 52 = s,$(s), B3 = SxF{s), (33) 

with n(s) = l//i(s), ^{s) = f2{s)/h{s), and F{s) = h{s)/h{s). 
From the last equation of ()32() we get 

u = (34) 

Substituting and into the system one can check that six equa- 
tions of (|32|) are satisfied identically and only the second equation leads to 

- 4:TrSx'^U{s)stt + 8TrstSxU{s)stx 

+ (4TrSx^U{s)p'p + Sx^F{sf - 4WU{s) + Sx^^sf) Sxx 



+ 47rs,.V. + 47rs,.Vn'(s) + s/^.(s)$'(s) + s/F(s)F'(s) = 0. (35) 
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Suppose that the pressure has the fohowing form 

^, , <t(sf + naf 2 

where G is an arbitrary function, then the equation 1)35(1 reduces to 

G'{s) 

sjstt - 2stSxStx + St^Sxx - "FFTT^ = 0- ("^6) 

n(s) 

It is possible to find two intermediate integrals ^ for the equation ()36|) 

St G'(s)t , , , St G'(s) 9 , , , , 

with (f) and ^p arbitrary functions. Eliminating St/sx from (|37() gives the im- 
plicit solution of ((36|l 

^ + '^(^)i + ^i' = V'(«). (38) 

In the special cases one can express s from (|38() and find the explicit solutions 
of the one-dimensional magnetohydro dynamics equations (|3'2() . 



4 Conclusion 

In this article we have developed the method of integrating system of first order 
partial differential equations in two independent variables. This method can 
be extended to the hyperbolic system of high order equations. It is important 
that corresponding invariants should exist for every family of characteristics 
of the system. 

When the system includes equations in n (n > 3) independent variables 
then we face a difficult task. As usual these systems have few invariants 
of characteristics. For example, the two-dimensional unsteady gas dynamics 
equations admit only one invariant of characteristics, namely the entropy. In 
three-dimensional case, the Ertel's integral is an additional invariant of char- 
acteristics. Note that in the case of the two-dimensional steady gas dynamics 
equations, Kaptsov have founded a first order invariant 

Jo — ) 

Sy 

where Ib is Bernoulli's integral 

Ib = + J -P,dp. 

Here, as usual, u and v are the components of velocity, p is the pressure, p is 
the density, and s is the entropy. 

There are problems that we have not yet studied. For example, it is in- 
teresting to look for new solutions of the one-dimensional gas dynamics equa- 
tions using the second order invariants (-^(4/5) ^(2)-'' ^^^^ interpretation to 
some founded solutions, and apply the Darboux method to other hyperbolic 
systems. To simplify calculations, we have implemented the package of ana- 
lytical computations which calculates characteristics and their invariants for 
the given system. 
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